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Abstract 

Bayes  estimation  of  the  reliability  function  for  the  inverse 
Gaussian  distribution  is  discussed.  For  the  case  that  the  mean  lifetime 
is  known,  Bayes  estimators  are  obtained  with  Jeffreys'  noninformative 
prior  and  with  the  natural  conjugate  prior  for  the  scale  parameter. 

In  the  case  that  both  parameters  are  unknown,  an  estimator  of  reliability 
is  suggested  which  is  based  on  the  Bayes  estimator  obtained  for  the  case 
that  the  mean  lifetime  is  known.  This  estimator  is  not  Bayes  but  compares 
favorably  with  the  maximum  likelihood  and  minimum  variance  unbiased 
estimators  as  indicated  by  computer  simulations. 

• 


-1- 


1  .  INTRODUCTION 

The  two-parameter  Inverse  Gaussian  distribution  with  probability  density 
function  (pdf)  in  the  form 

f(x|y,X)  -  (X/2itx^)^  exp  [-X(x-y)^/2y^x],  x  >  0,  p  >  0,  X  >  0,  (1.1) 

has  been  studied  in  the  reliability  or  life  testing  context  by  several  authors 

[6,  9,  10).  The  parameters  in  (1.1)  have  more  appealing  physical  interpretations 

for  life  testing  situations  than  other  parametric  forms  of  the  pdf.  The  mean  life 

for  the  lifetime  model  (1.1)  is  p,and  X  is  a  shape  parameter.  The  variance  is 
3 

p  /X  so  p  is  not  a  location  parameter  in  the  usual  sense.  Chhikara  and  Folks  [6] 
state  some  advantages  of  using  the  inverse  Gaussian  distribution  as  a  lifetime 
model  over  the  log-normal  distribution,  and  the  wide  variety  of  shapes  generated 
by  the  pdf  (1.1)  makes  it  a  competitor  to  other  lifetime  distributions.  In 
addition,  the  inverse  Gaussian  distribution  arises  as  the  first  passage  time  dis- 

i 

tribution  of  a  Brownian  motion  process  [7],  justifying  its  use  as  a  duration  time 
or  lifetime  model  on  a  physical  basis.  Several  results  have  also  been  obtained 
concerning  tests  for  drift  in  Brownian  motion  processes  (for  example  [5,  9,  13)). 

Tweedle  [15,  16)  and  Chhikara  and  Folks  [3,  4,  6],  among  others,  have  studied 
various  sampling  theory  inferences  concerning  (1.1).  Estimation  for  a  three-para¬ 
meter  Inverse  Gaussian  distribution  was  Investigated  recently  by  Padgett  and  Wei  [12]. 

The  cumulative  distribution  function  (cdf)  of  (1.1)  has  been  obtained  in 
closed  form  by  Shuster  [14]  and  Chhikara  and  Folks  [6],  and  the  survival  function 
or  reliability  is  given  in  the  form 

R(t|y,A)  -  ♦[(X/t)!*(l-t/|i)] 

-  exp(2X/y)  ♦[-  (X/t)-*  (1+t/p)  ],  t  >  0,  (1.2) 


where  $  denotes  the  cdf  of  the  standard  normal  distribution.  The  minimum 
variance  unbiased  estimator  of  R(t|p,X)  was  derived  by  Chhikara  and  Folks  C 6 J 
and  lower  confidence  bounds  for  (1.2)  were  given  in  [101.  However,  due  to  the 
complicated  nature  of  the  expression  (1.2),  other  inferences  concerning  the  re¬ 
liability  seem  to  be  difficult. 

Recently,  Banerjee  and  Bhattacharyya  [1]  presented  a  Bayesian  analysis  of 
the  inverse  Gaussian  distribution  in  a  different  parametric  form  than  (1.1)  (see 
Johnson  and  Kotz  [8]  for  other  forms  of  (1.1)).  As  was  stated  in  [11, 

Bayesian  inferences  concerning  reliability  are  extremely  difficult  and  require 
numerical  Integration  to  plot  the  posterior  pdf  or  to  determine  HPD  intervals 
for  reliability.  However,  Bayes  estimation  can  be  performed  in  come  cases,  and 
it  is  the  purpose  of  this  note  to  present  Bayes  estimators  of  R(t[p,X)  for  the 
case  that  the  mean  life  p  is  known,  which  is  reasonable  in  many  reliability  problems. 
Vague  priors  as  well  as  a  conjugate  family  of  prior  distributions  are  used.  In 
the  case  that  p  and  X  are  both  unknown,  an  estimator  of  (1.2)  is  propdsed  which, 
as  indicated  by  some  Monte  Carlo  simulation  results,  is  overall  as  good  as  the 
minimum  variance  unbiased  estimator  or  maximum  likelihood  estimator  given  in  [6] 
and  is  simpler  to  calculate  than  the  minimum  variance  unbiased  estimator.  The 
results  bear  a  remarkable  similarity  to  those  for  the  two-parameter  log-normal 
model  given  by  Padgett  and  Wei  [11]. 

2.  ESTIMATION  OF  RELIABILITY 

For  a  random  sample  x  -  (x^,  . ..,  x^)  from  the  inverse  Gaussian  distribution 

(1.1),  the  likelihood  function  is  given  by 

i  n  -3/2  X  x  n  -1 

i(X,p|x)  •  n  x.  1  exp  [-  ■=■  (==■  +  £  x.  )]  (2.1) 

i-1  1  i  u  i-1 

-1  n 

where  x  ■  n  £  x, .  The  Fisher  information  matrix  has  determinant 
i-1  1 
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1 1  (X  ,p  )  |  *  (w^X)  \  and  hence,  Jeffreys'  vague  prior  (Box  and  Tiao  T  2  J )  is 
p(p»X)  <*  (p  X)  2,  which  when  combined  with  the  likelihood  (2.1)  does  not  produce 
a  tractable  or  proper  posterior  distribution.  Also,  following  the  vague  prior 
idea  of  Box  and  Tiao  [2]  and  taking  p(p|X)“  constant  and  p(X)  «  X  \  mathemati¬ 
cally  intractable  posterior  distributions  for  estimating  R(t|p,X)  are  obtained. 

It  is  assumed  here  that  P ,  the  mean  life,  is  known,  and  Jeffreys’  noninformative 
prior  p (X )  «  X  ^  is  used  for  X.  In  addition,  the  gamma  family  of  distributions 
is  a  natural  conjugate  family  for  X,  and  Bayes  estimators  of  R(t|p,X)  for  this 
case  will  be  indicated. 

For  the  improper  prior  p(X)  *  X  \  the  posterior  distribution  of  X,  given  x, 
is  from  (2.1) 

p(X  |  x,p)  -  K  X  1  exp  [ - K-  E  (x.  -|i)2/x  1,  (2.2) 

2u  i-1 

where  the  constant  K  is  given  by 

.  n  j  ' 

K  -  r£)[-^  E  (x.  -U)7x.]. 

21i  1-1 

Hence,  p(X|jc,  p)  is  a  gamma  distribution  of  the  form 

p(X|x,p)  -  tr(a)Ba]  *  X°  1  exp(-X/B),  X  >  0, 

2  n  2 

with  a  ■  n/2  and  g  ■  2|i  /  1  t (x,  -  p)  /x.  ].  Then  with  respect  to  a  squared-error 

i-1  1  1 

*  2  n  2 

loss  function,  the  Bayes  estimator  of  X  ii  X.  ■  np  /  E  [(x.  -  p)  /x, which  is 

i-1  1  1 

the  same  as  the  mte  of  X  when  p  is  known.  For  R(t  |p,X),  with  respect  to  squared- 
error  loss,  the  Bayes  estimator  for  the  improper  prior  is  Rg(t)  -  EAtR(t  |  p,A)lxl. 
Thus,  from  (1.2),  for  each  t  >  0, 
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^(t)  -  Ea  t  $((A/t)’*  (1-t/y))  ) 

-  Eft[ exp(2A/|i)  ^(-(A/t) 2  (l+t/p))l.  (2.3) 

To  evaluate  the  first  expected  value  In  (2.3),  Lemma  1  of  Padgett  and  Wei 
fill  may  be  applied  with  c  In  that  lemma  equal  to  t  * ( 1  -  t/u)  .  Thua, 

EaI  ♦  ((A/t)^  (1+t/i.))]  -  PLT2a  <  (l-t/p)(a0/t)%),  (2.4) 

where  a  and  6  are  parameters  of  the  posterior  pdf  p(A|x,y)  defined  previously 
and  denotes  a  random  variable  having  Student's  t-distribution  with  v  degrees 
of  freedom.  The  second  expected  value  on  the  right-hand  side  of  (2.3)  is  evaluated 
similarly  after  absorbing  the  exponential  term  into  the  posterior  density  (2.2). 
Again,  by  Lemma  1  of  [11]  with  c  --t  ’  (1+t/p), 

EA[exp(2A/p)  ♦(  -(A/t)-4  (l+t/y))] 

-  (l-20/p)'a  P{T  <  -  (l+t/u)lo6|i/(t(p-23))]%).  •  (2.5) 

l  a 

Therefore,  the  Bayes  estimator  Rg(t)  is  given  by  (2.5)  subtracted  from  (2.4). 

This  estimate  may  be  easily  computed  since  it  involves  only  probabilities  for 
the  t  distribution. 

If  the  gamma  family  of  priors  with  parameters  y  and  6  in  the  form 
p(A)  «  A  exp(-A/6)  is  used  for  X,  then  the  same  kind  of  expected  values  are 

obtained  as  in  (2.3).  Applying  Lemma  1  of  [11]  again  yields  the  Bayes  estimator 
of  reliability  as 

RG(t)  -  PCT^  <  c1(y*/5*)J|] 

-  [«*/(6*-2/p)DT*  P<T2^  <  c2Cy*/(6*-2/w)]}, 


where  y*  -  y  +  n/2 


-5- 


<5*  =  *  *  +  < 2vi 2)  *  E  (y.j-pl^/x^,  c1  ■  t  *(l-t/p),  and  c^m -t  *(l+t/p). 

1  =  1 

These  results  bear  a  resemblance  to  those  of  the  log-normal  (or  normal) 
failure  model  obtained  in  [111  (see  also  [1J).  Also,  it  should  be  remarked 
that  for  the  noninforraative  prior  p(p|X)a  constant,  p(X)  “  X  * ,  estimates  of 
R(t|p,X)  may  be  obtained  by  numerical  integration,  but  a  closed-form  expression 
for  the  estimator  seems  extremely  difficult  to  obtain.  If  both  X  and  p  are 
unknown,  one  may  be  tempted  to  use  the  mi.e,  x  «*  p,  in  the  expressions  (2.4)  and 
(2.5)  to  obtain  an  estimate  of  reliability  R^(t).  The  effect  of  this  is  indicated 
in  the  next  section  by  some  computer  simulation  results. 

3.  MONTE  CARLO  SIMULATIONS 

Since  direct  comparisons  of  the  behavior  of  various  estimators  of  R(t|p,X) 
are  not  feasible  due  to  the  mathematical  complexity  of  the  estimators,  Monte  Carlo 
simulations  were  performed.  The  maximum  likelihood  (ML)  and  minimum  variance 
unbiased  (MVU)  estimators  when  p  is  known  were  compared  with  the  corresponding 
Bayes  estimator  (2.3).  For  several  values  of  t,  p,  and  X,  2000  samples  of  size 
n  (=  10,  20,  30)  were  generated  and  the  average  squared  errors  (ASE)  and  average 
estimated  reliability  (AER)  were  computed  for  each  estimator.  Similar  to  the 
results  in  (11],  the  Bayes  estimator  had  an  overall  smaller  mean  squared  error 
than  the  ML  and  MVU  estimators,  as  anticipated.  For  the  case  that  p  and  X  both 
were  unknown,  the  estimator  R^(t)  suggested  at  the  end  of  Section  2  (using  (2.3) 
with  p  replaced  by  x  -  p)  was  compared  with  the  ML  and  MVU  estimators  given  in  [6J 
in  the  same  kind  of  simulation  procedure.  Surprisingly,  this  estimator  performed 
as  well  as  the  MVU  estimator  in  the  sanse  of  average  value  and  did  not  have  a 
uniformly  larger  ASE  than  either  the  ML  or  MVU  estimator.  Some  of  the  results  of 
the  simulations  in  the  latter  case  are  given  in  Table  1 . 


Table  1.  Aver 


4 .  AN  EXAMPLE 


As  an  example,  the  estimator  R^(t)  as  well  as  the  ML  and  MVU  estimators 
were  used  to  estimate  reliability  for  several  values  of  t  from  the  n-46  repair 
time  observations  (in  hours)  for  an  airborne  communication  transceiver  (111!  and 
|6|).  Chlilkara  and  Folks  |6|  obtained  a  good  fir  to  this  data  by  the  inverse 
Gaussian  distribution  with  p  =  x  *  3.61  and  A  «  1.704.  The  estimates  of 
reliability  are  given  in  Table  2. 


Table  2.  Estimates  of  Reliability 


t 

1 

2 

3 

5 

10 

15 

s(t) 

0.6934 

0.4578 

0.3305 

0.1984 

0.0789 

0.0388 

MLE 

0.6986 

0.4607 

0.3325 

0.1996 

0.0791 

0.0386 

MVUE 

0.6951 

0.4618 

0.3368 

0.2057 

0.0829 

0.0396 

5.  CONCLUSION 

t 

For 

the  case  that  the  mean 

lifetime  p  in  the  inverse  Gaussian  model  is 

known,  the  posterior  distribution  of  A  is  easily  obtained  for  the  Jeffreys 
prior  and  the  natural  conjugate  prior  as  indicated  by  Banerjee  and  Bhattacharyya 
[1J.  For  this  case  the  Bayes  estimators  of  reliability  given  in  Section  2 
resemble  the  analogous  results  in  the  log-normal  (or  normal)  model.  If  both 
p  and  A  are  unknown,  the  Bayes  solution  for  reliability  in  a  compact  form  seems 
to  be  extremely  difficult,  at  least  for  the  parametric  form  (1.1).  It  also  seems 
to  be  even  more  difficult  to  obtain  a  Bayes  estimator  for  the  failure  rate 
function  or  mean  residual  life.  Hence,  an  estimator,  R^it),  of  reliability  was 
proposed  for  this  case  in  Section  2,  and  its  properties  were  Indicated  as  a  result 
of  computer  simulations.  For  other  Bayesian  inferences  on  reliability,  numerical 
integrations  must  be  performed  in  any  actual  application  to  obtain  the  posterior 
distribution  of  reliability. 
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